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^ A procedure  is  given  for  obtaining  an  exact  evaluation  of 
the  Rondo  integral.  This  is  a mathematical  function  used  to 
characterize  a type  of  electron  scattering;  it  arises  in  con- 
nection with  a variety  of  physical  processes.  Previously, 
there  was  no  exact  evaluation  of  the  integral,  and  various 
approximation  functions  were  used  to  represent  it.  Two  of 
these  approximations  are  compared  with  the  exact  answer. 
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1.  INTRODUCTION 


dilute  alloy  composed  of 


magnetic  Impurity  dissolved  in 


nonmagnetic  host  may  exhibit  anomalous  resistivity  at  low  temperature. 
Kondo1  explained  the  phenomenon  by  taking  into  account  the  spin-spin 
interaction  between  conduction  electrons  and  immobile  ions.  He  showed 
that  this  led — in  third-order  perturbation  theory — to  a resistivity 
having  logarithmic  dependence  on  temperature.  The  same  type  of 
electron-ion  coupling  has  been  used  to  explain  other  phenomena.  For 
example,  Kondo  scattering  has  been  invoked  to  explain  anomalous 
tunneling  in  thin-film  structures 2 and  Schottky-barrier  junctions,3  and 
to  explain  negative  magnetoresistance  in  semiconductor  materials.4 

Until  now,  no  exact  evaluation  of  the  Kondo  integral  has  been 
reported.  Applebaum5  performed  a numerical  integration  and  suggested  an 
approximation  function  to  represent  his  results.  Wolf  and  Losee6  also 
did  a numerical  evaluation,  but  they  suggested  a different  approximation 
function.  The  work  reported  here  gives  an  exact  evaluation  of  the 
integral  in  terms  of  analytic  functions. 


2.  FORMULATION 


Let  the  Kondo  integral  be  defined  as 


where  f(x)  is  the  Fermi  function.  The  function  g(x,c)  is  defined  as 


g(x,c) 


where  £?  denotes  principal  value  of  the  integral,  and  c is  a cutoff 
parameter . Equation  (2)  can  be  integrated  by  parts.  Proper  care  must 


lJ.  Kondo,  Prog.  Theoret.  Phys.  (Kyoto)  32_  (1964),  37. 

2J.  Appelbaum,  Phys.  Rev.  Lett.  17  (1966),  91. 

3P.  H.  Anderson,  Phys.  Rev.  Lett.  17~  (1966),  95. 

**R.  P.  Khosla  and  J.  R.  Fischer,  Phys.  Rev.  2B_  (1970),  4084 

5J.  Appelbaum,  Phys.  Rev.  154  (1967),  633. 

6ff.  L.  Holf  and  D.  L.  Losee,  Phys.  Rev.  B2  (1970) , 3660. 


be  taken  with  the  singularity 
£n|t  - x|  is  the  integral  of  (t 
Therefore,  the  result  is 


at  t = x,  and  it  is  noted  that 
x)-1  dt  for  both  t > x and  t < x. 


g(x,c)  - -[f (-c)-*s]*.n|c  + x|  + [f (c)-y in|c 
- /Cf'  (t) in |t  - x|dt  . 


The  next  step  is  to  substitute  equation  (3)  into  equation  (1)  and 
perform  the  integration.  For  this  it  is  useful  to  have  the  result 
proved  in  appendix  A,  namely 

J (a)  - -y"f(t)in|t  - a|dt  - Intfir)  + Re  ¥<»j  + ia/2ir)  . (4) 


In  equation  (4)  Re  denotes  the  real  part  of  the  complex  function, 
and  ¥ is  the  digasua  function.7  Using  equation  (4),  one  can  substitute 
equation  (3)  into  equation  (1)  and  obtain  the  following: 


F(y,c)  ■ -[f (-c)-J,]{in2x  + Re  ¥[«*  - i(c  + y)/2ir]} 


+ lf(c)-H]Un2ir  + Re  ¥[>,  + i(c  - y)/2x]> 


i(t  + y)/2n] }dt 


In  the  problems  of  physical  interest,  such  as  tunneling  and  ancsialous 
resistivity  in  alloys,  the  parameter  c has  a magnitude  of  the  order  of 
100.  In  such  cases,  it  is  an  excellent  approximation  to  let  f(-c>-*-l, 
f(c)->-0,  and  use  the  large-argument  approximation  for  the  digamna 
function,  i.e.. 


With  these  approximations,  equation  (5)  reduces  to  the  form 


i(t  + y)/2x]dt 


M.  Abramowitz  and  I.  Stegun,  Handbook  of  Mathematical  Functions 
Government  Printing  Office , Waahington,  D.C.  (1964),  2S8. 


Appendices  B and  C show  how  to  perform  the  integration  in  equation  (7) . 
For  c -►  °»,  the  integration  can  be  performed  by  using  the  theory  of 
residues  (shown  in  appendix  B) . Then  a correction  term  is  derived  for 
finite  c (shown  in  appendix  C) . For  any  problem  of  physical  interest, 
the  correction  term  is  very  small  because  the  major  contributions  to  the 
integral  in  equation  (7)  come  from  those  regions  near  t ■ 0 and  t ■ -y . 
The  integrand  decreases  rapidly  for  large  t,  and,  if  c is  large , the 
integral  has  almost  the  same  value  for  finite  c as  for  c -*-<*>.  Using 
the  results  developed  in  appendices  B and  C,  one  can  write  the  exact 
evaluation  of  the  Kondo  integral  as 


F(y,c)  - -hln \ c2  - y2!  + in(2ir)  - 1 + Re  ¥(1  + iy/2ir) 


+ yr=-  Re  ¥(1  + iy/2n)  + K(y,c)  . (8) 
dy 

As  before.  Re  denotes  the  real  part  of  the  complex  variable,  and  \|/  is 
the  digamma  function.  The  term  K(y,c)  in  equation  (8)  is  the  correction 
term  for  finite  c that  was  developed  in  appendix  C. 


A discussion  of  the  Kondo  function  follows  in  section  3 of  this 
report.  First,  however,  a remark  may  be  made  concerning  the  function 
g(x,c)  given  in  equation  (3).  Certainly  the  principal  dependence  of  g 
on  c is  embodied  in  the  first  two  terms  of  equation  (3) ; little  error  is 
incurred  by  letting  the  limits  on  the  integral  become  infinite. 
Therefore  one  can  use  equation  (4)  to  write 


g(x,c)  » - l n(c2  - x2)  + in(2ir)  + Re  ¥ (**  + ix/2ir)  . 


This  agrees  with  the  result  given  by  Bloomfield  and  Hamann,8  (their 
equation  A-8) , but  they  "derived"  the  digamma  function  by  deducing  a 
recurrence  relation,  and  then  arguing  that  the  recurrence  relation 
defined  the  function  to  within  an  additive  constant.  Actually,  the 
recurrence  relation  of  the  digamma  function  defines  it  only  to  within  an 
additive  function,  since,  for  example,  i|>  (z)  + A sin  (irz)  obeys  the  same 
recurrence  relation  as  'Hz). 


®p.  Bloomfield  and  D.  R.  Hamann,  Phya.  Rev.  164  (1967),  856. 
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3.  APPLICATION  OP  RESULTS 


Table  I gives  numerical  values  of  the  function 
H(x)  - £n(2w)  - 1 + Re  f ft  + ix/2ir)  +x^Reftt  + ix/2ir)  . 


U0> 


TABLE  I.  NUMERICAL  VALUES  FOR  THE  FUNCTION 
H(y)=tn(2ir)-1+Re  f (l+iy/2ir)+yd/dy 
Re  T(l+iy/2ir) 


y 

Re  f 

d/dy ( Re  Y) 

H (y ) 

0.0 

-0.5772 

0.0000 

0.2607 

0.1 

-0.5769 

0.0382 

0.2616 

0.2 

-0.5760 

0.0764 

0.2643 

0.3 

-0.5745 

0.1143 

0.2689 

0.4 

-0.5724 

0.1520 

0.2752 

0.5 

-0. 5696 

0.1892 

0.2833 

0.6 

-0.5663 

0.2260 

0.2931 

0.7 

-0.5625 

0.2622 

0.3046 

0.8 

-0.5580 

0.2977 

0.  3178 

0.9 

-0.5530 

0.3325 

0. 3325 

1.0 

-0.5474 

0.3665 

0.3488 

1.1 

-0.5413 

0.3996 

0.3665 

1.2 

-0.5347 

0.4317 

0.3856 

1.3 

-0.5276 

0.4628 

0.4061 

1.4 

-0.5200 

0.4929 

0.4277 

1.5 

-0.5119 

0.5219 

0.4506 

1.6 

-0. 5034 

0.5497 

0.4745 

1.7 

-0.4944 

0.5763 

0.4994 

1.8 

-0.4850 

0.6017 

0.5252 

1.9 

-0.4753 

0.6259 

0.5519 

2.0 

-0.4651 

0.6489 

0.5793 

2.  1 

-0.4546 

0.6706 

0.6074 

2.2 

-0.4438 

0.6911 

0.6361 

2.3 

-0.4326 

0.7103 

0.6653 

2.4 

-0.4212 

0.7283 

0.6949 

2.5 

-0.4094 

0. 7451 

0.7249 

2.6 

-0. 3975 

0.7606 

0.7552 

2.7 

-0. 3852 

0.7750 

0.7857 

2.8 

-0.3728 

0.7882 

0.8163 

2.9 

-0. 3602 

0.8003 

0.8471 

3.0 

-0. 3473 

0.8113 

0. 8779 

3.1 

-0.3343 

0.8212 

0.9087 

3.2 

-0.3212 

0.8301 

0.9394 

3.3 

-0.3079 

0.8380 

0.9701 

3.4 

-0.2945 

0.8449 

1.0006 

3.5 

-0.2810 

0.8509 

1.0  30  8 

3.6 

-0.2674 

0.8560 

1.0609 

3.7 

-0.2538 

0.8603 

1.0907 

3.8 

-0.2401 

0.8638 

1.1202 

3.9 

-0.2263 

0. 8665 

1.1494 

4.0 

-0.2125 

0.8685 

1.1783 

r;v 


■ 
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This  is  the  Kondo  function  without  the  terms  that  depend  on  the 
parameter  c.  Values  in  table  I were  calculated  from  a series 
representation  of  the  diyamma  function,6  i.e., 

n 

Re  ¥(«*  + iy)  « -y  + Z y2/n(n2  + y2)  . (11) 

n=l 

The  sumations  were  continued  until  the  last  term  was  less  than 
10~7  times  the  sum  of  the  preceding  terms;  the  same  rule  was  used  for 
calculating  the  derivative  of  equation  (11 ) . 

Figure  1 compares  the  correct  Kondo  function  and  the  approximation 
used  by  Appelbaum. 5 His  approximation  function  was 

F approx  = -iin[(y  + l)/c]  . (12) 


•-  v 


APPROXIMATE 


Figure  1.  Kondo  function  for  £0  = 10  meV  and  kT  = 0.1  meV, 
corresponding  to  c ■ EqAT  = 100.  The  dashed 
curve  is  the  approximation  used  by  Appelbaum. 


5J.  Appelbaum,  Phys.  Rev.  154  (1967),  633. 

8P.  Bloomfield  and  D.  R.  Hamann,  Phys.  Rev.  164  (1967),  856. 


This  function  has  finite  slope  at  y = 0,  whereas  the  correct  Kondo 
function  has  zero  slope.  The  value  of  the  parameter  used  in  figure  1 

was  c ■ 100,  corresponding  to  Appelbaum ’ s cutoff  energy  E0  = 10  meV  at  a 
temperature  kT  * 0. 1 meV ; thus  c = E /kT  = 100.  The  Kondo  function  has 
smaller  slope  than  the  approximation  function,  especially  as  y -*•  0. 
This  difference  could  be  important.  For  example,  in  the  theory  of 
magnetoresistance , y is  the  ratio  of  magnetic  energy,  yB,  to  thermal 
energy,  kT.  For  a magnetic  moment  of  one  Bohr  magneton,  in  a field  of 
one  kilogauss,  at  the  temperature  of  liquid  helium,  y = yB/2irkT  = 10" 3 . 
This  is  small  and,  furthermore,  one  uses  the  derivative  of  the  Kondo 
function  to  analyze  magnetoresistance  data.  Misleading  results  would  be 
obtained  if  one  used  the  approximation  function  instead  of  the  correct 
function  in  this  case. 

Figure  2 compares  the  correct  Kondo  function  with  the  approximation 
function  used  by  Wolf  and  Loses.6  Their  approximation  function  was 

F approx  - -in[(y2  + a2)1*/^]  . (13), 


Figure  2.  Kondo  function  for  EQ  - i meV  and  T * 1.25  K,  corresponding 
to  c * Ejj/kT  * 18.55.  The  dashed  curve  is  the  approximation 
function  of  Wolf  and  Losee  with  a ■ 2.12. 


'1 


t| 


This  function  has  the  correct  slope  at  y **  0,  but  the  wrong  magnitude  at 
y - 0 unless  a is  chosen  properly.  Wolf  and  Losee  used  the 
value  a ■ 2.12,  which  causes  F(y)  to  be  too  large  at  y * 0.  Also,  the 
correct  Kondo  function  has  larger  slope  than  the  approximation  function 
for  all  y > 0. 

Shen  and  Rowell3  stated  that  the  result  of  a numerical  integration 
of  the  integral 

00 

ftne  =/ / f ' (t)f  • (x)  in  1 1 - x | dx  dt  (14) 

—00 

was  3 = 1.35.  The  exact  result,  obtained  from  appendix  B,  is 

inB  = in(2ir)  - 1 - y . (15) 


This  gives  the  correct  value,  3 = 1.2978. 

Problems  of  physical  interest  frequently  require  an  evaluation  of 
the  Kondo  function  at  values  of  the  argument  small  compared  with  unity. 
In  such  cases,  it  is  convenient  to  use  an  approximation  of  the  digamma 
function  derived  from  equation  (11),  i.e.. 

Re  y (*a  + iy)  = - y + y2  C ( 3)  , (16) 

where  ^(3)  * 1.202  is  the  zeta  polynomial.  Then  the  Kondo  function  can 
be  written 

f (y,c)  = -inc  + a + by2  (17) 

where  a = £n(2ir)  - 1 - Y = 0.2607,  and  b = 3C(3)/(2tt) 2 = 0.0776. 


4.  SUMMARY 

The  Kondo  function,  defined  in  equations  (1)  and  (2) , was  evaluated 
by  contour  integration  and  shown  to  be  expressible  in  the  form  given  in 
equation  (8) . This  differs  slightly  from  the  approximation  functions 
used  by  various  authors.  For  example,  the  correct  function  has  zero 
slope  at  the  origin,  whereas  the  function  suggested  by  Appelbaum  had 
finite  slope.  The  function  used  by  Wolf  and  Losee  had  correct  slope  at 
the  origin,  but  the  magnitude  was  incorrect.  In  general,  the  correct 
function  has  larger  slope  than  the  Wolf  and  Losee  function  and  smaller 
slope  than  the  Appelbaum  function,  for  y •>  0.  A definite  integral 
mentioned  by  Shen  and  Rowell  was  evaluated  exactly  and  found  to  have  a 
smaller  value  than  the  approximate  value  they  quoted. 


*L.  Y.  L.  Schen  and  J.  M.  Rowell,  Phys.  Rev.  165  (1968),  566. 
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APPENDIX  A 


The  function  f ' (z)  has  second-order  poles  at  z = iir  (2n  ^1) , where 
n * 0,  1,  2,  . . . , and  the  residue  of  the  integrand  at  the  n pole  is 


R - [a  - iir  (2n  + l)]"1  . 
n 


(A -4) 


The  contours  encloses  N poles,  and  the  procedure  is  to  take  the  limit  as 
N -*■  00 , after  the  integral  has  been  evaluated  for  finite  N. 


The  integration  consists  of  four  parts:  integration  along  the  real 
axis,  Ij,  integration  along  the  top  contour  parallel  to  the  real  axis, 
l2>  and  integration  along  the  two  vertical  legs.  The  latter  contribute 
nothing  because  f ' (z)  ■*  0 as  the  real  part  of  z -*■  ±°°.  Therefore, 
according  to  the  theory  of  residues,  the  result  can  be  written  as 

N-l 

Ii  + I2  * 2iri£  R . (A -5) 


n=0 


Along  the  real  axis,  the  phase  angle  is  <J>  ■ n for  u < a,  and  4>  = 0 
for  u y a.  Substitution  of  equation  (A-3)  into  equation  (A-l)  therefore 
leads  immediately  to  the  result 


II  = ,J  (a)  - iir[f(a)-l]  , 
where  J(a)  is  the  integral  along  the  real  axis,  i.e. 

J(a)  - -J  f'(u)£n|u  - a|du  . 

— CO 

Along  the  upper  contour,  the  logarithm  can  be  written 
ln(u  + i2Nir  - a)  * £n(i2Nir)  + £n[l  + (u  - a)/i2irN] 


(A -6) 


(A -7) 


(A-8) 


The  second  term  in  equation  (A-8)  is  less  than  (u  - a)/i2Nn,  which 
approaches  zero  as  N -►  ®.  Therefore,  the  result  of  integration  along 
the  upper  path  is 


12  - -in(i2Nir)  + 0(1/N)  . 
Using  these  results,  one  has 


(A-9) 


N 


J (a)  - in  [f  (a)  - 1]  - in  (2iNir ) - 2ni  £ [a  - in  (2n  + 1)]  . (A-10) 

n=0 
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, 


Let  K « *s  + i (a/2* ) 
can  be  written  as 


Then  the  term  on  the  right  side  of  equation  (A-10) 


Using  the  recurrence  relation 


J<€)  - iir[f(a)-y  + tn2ir  + Y(C) 


Note  that  the  imaginary  part  of  Yt’s+iy)  is  (ir/2)  tanh  Try,  which  means 
that  Xmi//(£)  exactly  cancels  the  first  term  in  equation  (A-15).  The 
final  result  is 


J(a)  = tn2ir  + Re  iM*!  + (ia/2ir)] 


1M.  Abramowi tz  and  I.  Stegun,  Handbook  of  Mathematical 

Functions,  U.S.  Government  Printing  Office,  Washington,  D.C.  (1964 J , 
259. 
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[Re  t(l  + ix/2ir) ) 


H(x)  - -1  + Re  # (1  + ix/2ir)  + x 


function  is1 


1 - eotp(-u) 


The  procedure  is  to  substitute  equation  (B-3)  into  equation  (B-l) , 
letting  s ■ **  + i(t  + x)//2ir.  one  can  perform  two  trivial  integrations 
immediately.  Then,  having  interchanged  the  order  of  integration  in  the 
third  term,  one  obtains 


f’(t)exp(-i  t/2w)dt  dw 


The  t-integral  in  equation  (B-4)  can  be  done  by  contour  integration 
Consider  the  function 


defined  by  the  contour  shown  in  figure  B-l.  The  denominator  of 
equation  (B-5)  has  second-order  poles  at  r ■ iri(2n  + 1) , 
n ■ 0,1,2  . . .,  and  the  oontour  encloses  one  of  these  poles.  The 
residue  at  z ■ iri  is  R * i a exp(ira).  Let  the  integral  along  the  real 
axis  be  Xi,  and  the  integral  along  the  upper  leg,  parallel  to  the  real 
axis,  be  X2.  X2  is  the  same  as  Xi,  except  that  the  variable  in  I2  is 
s • u ♦ 2 ni,  instead  of  z • u,  as  it  is  in  X}.  This  introduces  only  a 
multiplicative  factor,  so  that 


1M.  Abramcsritz  and,  I.  Stegun,  Handbook  of  Mathematical 
Functions,  V.S.  Government  Printing  Office , Washington , D.C.  (1964), 
259. 


APPENDIX  B 


/v 


Figure  B-l.  Contour  in  the  complex  plane  used  in 

evaluating  the  integral  of  equation  (B-l) . 


The  integrals  along  the  vertical  legs  do  not  contribute  anything  because 
the  integrand  in  equation  (B-5)  approaches  zero  as  u -*■  ±®.  Therefore, 
according  to  residue  theory,  the  result  of  the  contour  integration  is 

II  II  - exp  (2ira)  ] - T2rra  exp(ira)  , (B-7) 


which  becomes 


II  (a) 


tra 

sinh(ira)  ' 


CB-8) 


Except  for  a minus  sign,  the  t-integral  in  equation  (B-4)  is  the  same  as 
Il(a),  with  a - x/2ir.  Therefore,  equation  (B-8)  can  be  substituted  into 
equation  (B-4 ) , and  the  result  written  as 


H(x) 


exp  (-u)  du 
1 - exp(-u) 


+ 


Re 


uexp[-u(l  + ix/2ir))du 
[1  - exp  (-u) ] 2 


(B-9) 


(B-10) 


t (u)  - u exp  (1  - ux/2it) 


g (u)  - -exp(-u)/[l  - exp(-u)] 


exp(-u)du 
1 - exp(-u) 


By  using  (B-3),  one  can  rewrite  this  as 


H(x)  ■ -1  + Re  i(i{1  + ix/2ir) 


[Re  <|>(1  + ix/2ir) ] 


This  is  the  result  claimed 


APPENDIX  C . — CORRECTION  TERM  FOR  FINITE  CUTOFF 


This  appendix  will  show  that  the  error  incurred  by  letting  c -*■  ® In 
equation  (5)  of  the  body  of  this  report  is  acceptably  nail.  The  amount 
added  to  the  integral  when  the  upper  lladt  approaches  infinity  is 


f'(t)  Re  ¥[Jj  + iCt  + x)/2ir]  dt 


K(x,c)  < f exp(-t)ln[(t  + x)/2ir]dt 


K(x,c)  < exp(-c)ln[(C  + x)/2ir]  + f exp(-t)(t  + x)-1dt 


(C-6) 


Therefore,  the  final  result  can  be 


K(x,c)  < expl-c)  {1  ± lnl(t  + x)/2tt]> 


Sample  values  of  the  correction  term  are  K(5,10)  > 4.6  * 10 
K(10,10)  < 5.3  x io~9.  These  are  small  compared  with  the  values 
corresponding  Integrals. 


‘if.  Abranowitz  and  I.  Stegun,  Handbook  of  Mathematical 
Functions,  Dover  Publications,  Inc.  (1965) , 25 9. 
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